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Abstract 

We discuss two dedicated observables in exclusive B K*i'^i~ decay that can be 
used to study effects of physics beyond the standard model, namely the forward- 
backward asymmetry in the lepton spectrum and the isospin-asymmetry between 
decays of charged and neutral B mesons. We consider the region of large recoil- 
energy (i.e. small invariant mass of the lepton pair), and employ the QCD factoriza- 
tion approach to exclusive B meson decays. Sub- leading effects in the heavy quark 
mass expansion have been taken into account for the calculation of the isospin- 
asymmetry. We give predictions for decay asymmetries in the standard model, 
and its supersymmetric extension with minimal flavor violation, using parameter 
values allowed by current experimental constraints on — > Xg^ decay. 



1 Introduction 



B physics has become one of the most active fields in the search for new interactions 
associated to physics beyond the standard model (SM). In particular, the fiavour sector 
of the new theory should have distinctive effects on exclusive and inclusive rare B decays 
which are suppressed in the SM. At present, one of the most important constraints on 
new physics from rare B decays arises from the inclusive decay B X^'y. 

The inclusive decay processes B Xs'y and also B Xal^l~ have been studied 
in great detail in the literature in the SM and MSSM 0-||16[. Concerning B — >• 

we have at our disposal accurate determinations both from experimental [|l^ and 
theoretical point of view [0-0]. Within the present uncertainties B —>■ Xsj is consistent 
with the SM prediction, which translates into a quite narrow allowed region for the 
Wilson coefficient C^^ in the weak effective Hamiltonian for b ^ s'j (also b — > 577 p8[| ) 
and b — > si~^i~ transitions. The inclusive B Xsl^l~ decay has not been measured 
yet; for the theoretical literature we refer the reader to [§|,|T^-[0. 

In this paper we will focus on the exclusive decay mode B K*l^l^ . In general, 
exclusive decays are much easier accessible in experiment (first experimental evidence 
for these decay mode has been presented in ^), but their theoretical treatment is 
usually more difficult because of non-perturbative hadronic binding effects. However, 
in the limit where the invariant mass of the lepton pair is small, the process admits a 
systematic theoretical description using QCD factorization in the heavy quark limit 



As stressed in |2^, , a distinct observable in this decay is the forward-backward 

(FB) asymmetry in the lepton spectrum: In particular the FB asymmetry-zero can be 
calculated with rather high accuracy since most of the hadronic uncertainties drop out. 
This leads to an almost model-independent relation between the Wilson coefficients C^^ 
and Cg. 

We further stress that the sensitivity to new physics in exclusive modes may be 
different from that of their inclusive counterparts if the spectator quark is involved 
in the short- distance dynamics. This is the case for the isospin asymmetry between 
B K*i^i~ decays with charged or neutral B mesons. Kagan and Neubert have 
shown that annihilation topologies provide the main source for the isospin asymmetry 
in B —>■ K*'y decays, and concluded that the sign and magnitude of the isospin asym- 
metry is predominantly determined by the Wilson coefficients C5 and Cq of four-quark 
penguin operators and their relative sign with Cf^ . The SM prediction is consistent with 
experiment, within still large uncertainties. In this article we will extent the analysis of 
p6| to the decay B — > K*i~^i~ and calculate the isospin-asymmetry for small invariant 
lepton-pair masses (g^ below 7 GeV^). 

The paper is organized in the following way. In the next section we will shortly review 
the FB asymmetry for the decay B K*l'^l" within the QCD factorization approach. 
Then we turn to the calculation of the isospin asymmetry in 5 — >■ K*l~^l~ , and give a 
numerical estimate within the SM. In section 3 we explore new physics contributions to 
both observables, the FB and the isospin-asymmetry, in the MSSM with minimal flavor 
violation (MFV). We apply the constraints from B —>■ Xgj using the NLO calculation 
in [|T^, |TI|. Section 4 is devoted to a brief discussion of possible effects beyond MFV, 
before we conclude with a summary in section 5. 
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2 B ^ K*i^i decay asymmetries in the SM 

The decay B K*£^£~ is induced by a set of operators Oi appearing in the weak 
effective Hamiltonian for b ^ s transitions [HI 



^ r 2 10 ^ 

^eff = ^ |E(^« or + A, a o'r) - E a o,| 



(1) 



Here, the combinations of CKM matrix elements are abbreviated as Xq = V*^ Vgh, and 
obey the unitarity relation A„ + Ac = — At. The Wilson coefficients contain the short- 
distance physics related to integrating out heavy particles (like weak gauge bosons and 
the top quark) and hard gluons. In extensions of the SM new (heavy) particles modify 
the SM values of Wilson coefficients and/or lead to new operators in addition to (|l|). In 
this paper we will restrict ourselves to the case of SM operators only. (We will use 



the same conventions as in in particular, the gauge vertices read +ieeg and +igsT , 
and the SM values for Cf and Cf are negative.) 

A systematic calculation of QCD effects in exclusive B decays requires the factor- 
ization of long- and short-distance physics distinguished by the heavy quark mass rrib in 
the hadronic matrix elements of Hes- The theoretical framework has been introduced 
for non-leptonic B decays by Beneke/Buchalla/Neubert/Sachrajda (BBNS) and 



has been extended to heavy-to-light form factors and radiative B decays in [E^, |EB 



| 3T| . In the BBNS approach, starting from the heavy quark limit, rrib oo, radiative 
corrections to "naive" factorization of hadronic amplitudes arising from hard gluon ex- 
change can be calculated in terms of perturbative hard-scattering kernels and universal 
(non-perturbative) hadronic quantities (decay constants, form factors, and light-cone 
distribution amplitudes). 

In order to constrain new-physics scenarios in exclusive decays with reasonable ac- 
curacy, one has to find observables that are sensitive to the short-distance physics and 
rather independent to hadronic input parameters like form factor values, moments of 
light-cone wave functions etc. In this paper we will focus on two asymmetries in the 
lepton spectrum, namely the forward-backward asymmetry 



~ ' ' d{cosO) / d[cosP^ 



dq^ dV/dq^ \Jo dq^d cos 6 J-i dq'^d cos 6 

and the [CP averaged) isospin asymmetry 

dAi _ rfr[fiO ^ K*H+i-]/dq^ - dT{B^ K*^i+i-]/dq^ 



dq^ dr[BO K*H+i-]/dq^ + dr[B^ K*^i+i-]/dq'^ 



(2) 



(3) 



Because of the sensitivity to relative signs and orders of magnitude between different 
Wilson coefficients, these asymmetries are particularly useful to constrain new physics 
scenarios. (|^ entertains another possibility to study physics beyond the SM in this B 



decay mode, by comparing helicity amplitudes for B K*£^£ and B — > pi^u.] 
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2.1 Forward-backward asymmetry 



We will follow the results and notation of [£2[ where the QCD corrections to the decay 
B — > K*£~^i~ have been calculated to leading power in the inverse heavy quark mass 
or light-quark energy, respectively, and to NLO in the strong coupling constant. The 
approach is restricted to the situation where the energy of K* is large, i.e. the invariant 
lepton-pair mass is small. For practical purposes this means that is below charm 
threshold, < 4m^ ~ 7 GeV^. The decay amplitude for B —>■ K*i~^i~ is conveniently 
expressed in terms of the quantities Cg{q'^) and Cg(g^), which are generalizations of the 
effective Wilson coefficient Cf^ used in the "naive" factorization approach. The lowest- 
order (naively factorizing) expressions read 



•-9 



2^ , 2mbMB^,ff 



2m5 



Cf . 



B 



(4) 



(The function Y{s) can be found, for instance, in [^, ^). The differential decay rate 



then reads 22 



dq^dcos 6 



128713 M^X^q.m^,) y 



+ (1 - cos^ 



-^o,9j^^U{q 




\Cl{q')? + ClA^\{q^f) 



2\2 



Mc^iq'))C- 
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Here 



Hq' 



m 



K* 



[(^ Ml) 


2ml, 








Ml 







1/2 



(5) 



(6) 



is the usual kinematic factor, m^ps is the potential-subtracted heavy quark mass [Q, 
E ^ {Ml — q^)/2MB is the K* energy, and 6 refers to the angle between the positively 
charged lepton and the B meson in the center-of-mass frame of the lepton pair. Further- 
more ^±{q^) and ^||(g^) are the universal soft heavy-to-light form factors, arising in the 
heavy-quark and large-recoil-energy limit |3^, ^ (for a critical examination whether the 
form factors themselves can be calculated perturbatively, see and references therein.) 
The function A||(g^) is given in |^ . 

From (^ it is obvious that the FB asymmetry dA-pB/dq^ in (H) is proportional to 
Re {C^{q^)^, and therefore it vanishes if Re (C^(go)) = 0- leading order this gives a 
relation between the short- distance Wilson coefficients Cg and Cj^ which is completely 



independent of hadronic uncertainties ||2^, 

LO: C^ + Re{Y{ql)) 



2MBmb 

ql 



■eff 



(7) 
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Figure 1: The SM forward-backward asymmetry dA-pB,/dq^ for the decay B~ —>■ K*~i~^i~ 
as a function of the invariant lepton-pair mass at next-to-leading order (sohd center 
hne) and leading order (dashed). The band reflects all theoretical uncertainties from 



input parameters and scale dependence combined. Figure taken from [22 



Therefore the measurement of the asymmetry zero provides a crucial test of the SM 
and new physics. (Phenomenological consequences for the FB asymmetry from physics 
beyond the SM have already been discussed in using "naive" factorization.) In Fig. |I| 
we plot the results for the FB asymmetry obtained in ||22[ for SM Wilson coefficients. 



A conservative estimate for the location of the asymmetry zero in the SM (taking into 
account the possibility of additional model-dependence from the form factors arising at 
higher order in the Ah/ms expansion) was obtained as 



(4.2 ± 0.6) GeV^ 



(8) 



2.2 Isospin asymmetry 

The isospin-asymmetry arises from "non-factorizable" graphs where a photon is radiated 
from the spectator quark in annihilation or spectator-scattering diagrams, see Fig. ^ 
The corresponding contribution to the partial decay width depends on the charge of the 
spectator quark, and therefore it is different for charged and neutral B meson decays. 
The isospin asymmetry vanishes in naive factorization. Non-factorizable annihilation 
and spectator-scattering topologies induce only a tiny asymmetry when working in the 



heavy quark limit. However, it has been observed in l^g that for B — > K*'y decays 
isospin-breaking effects that formally are of sub-leading order in the Ah/rriB expansion 
numerically may be large. Although part of the calculation requires to model new non- 
perturbative effects, the dominant isospin-breaking effect stems from processes which are 
dominated by short-distance effects and are thus perturbatively calculable in the QCD 
factorization approach. 



In the following, we extend the analysis in [EH] to the B — ^ K*i^i case, including 



also longitudinal photon polarizations appearing for q ^ 0. For that purpose we consider 



4 



B 




-K* 



B 




K* 



(b) 



B 



(c) 



B 




K* 



B 



K* 



B 




K* 



-K* 



Figure 2: Contributions to the isospin asymmetry dAj/dq"^: (a) Annihilation topolo- 
gies with operators Oi_6, (b) Hard spectator interaction involving the gluonic penguin 
operator Og, (c) Hard spectator interaction involving the operators Oi_6- 



the exclusive quantities Cg and Cg introduced in |Q, and define 



where the leading-order expressions Cg^^^iq'^) and C9°^"(g^) are given in and the 
functions hq{q^) and h^\{q^) parameterize non-factorizing effects from photon radiation 
off the spectator quark line with flavor q. Non-factorizable effects that are independent 
of the spectator quark will drop out of the isospin asymmetry to the considered order of 
the as and l/mi, expansion. (Note that there is a formally LO annihilation contribution 
to Cg which does depend on the charge Cq of the spectator quark and will therefore 
appear in the function Because this contribution is numerically small we will 

treat it as a sub-leading effect) In terms of these new parameters the (CP-averaged) 
isospin asymmetry (§) can be written as 



(9) 



dAi 
dq^ 



Re (6, 
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|/-(0)- 
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ef(g') Re(6! 



q'^mj^. 



aiq' 



\c 



(0)11,2 
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|(-{0)±|2 



-1 
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IC, 
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+ \c 



10 1 



(10) 



where we neglected terms of order 6^. In the limit 



the photon-pole in Cg^^ ^ will 
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dominate, Cg'^' 
Aj[B ^ K*^] 



l/q^i see @). In this hmit the isospin-asymmetry in ([lOD reduces to 



Re 



(0) ~ (0) 5 which is the quantity discussed in For larger 



values of the isospin asymmetry will be dominated by the longitudinal polarization. 



Generalizing the definitions in p6|, we further define the functions ^-^d Kl via 



bliq') 



24 vr^ fs eg f /a-* 



+ 



Here 



duo 



6\B,+ {q^)mB l-q^/m\ 
-^Kl{q')] ■ 



(11) 



(12) 



/o UJ — q^ /ruB — ie 
are moments of the B meson light-cone distribution amplitudes as defined for instance 



m 



36, 27 . The function K-l receives an annihilation contribution to leading order in the 



^h/f^B expansion with 



(13) 



Here the Wilson coefficients Ci refer to the basis defined in |^ (at leading order it 
coincides with the basis used in 0). The functions in (O) only appear at sub- 
leading order of the A/i/m^ expansion. Nevertheless, numerically the leading order 
effect in lA{q^) was found to be rather small ||2^, whereas the formally sub- leading effects 



entering &g (0) were found to be sizeable [^. In particular, the contributions from the 
scalar penguin operators O5 and Oq play an important role, despite the fact that the 
corresponding Wilson coefficients are rather small. In the following we will extend the 
calculation of Kagan/Neubert in and study how these effects change with increasing 
photon virtuality q^. Note that in this section we will only consider those corrections 
that are sensitive to the charge of the spectator quark, and neglect all sub-leading effects 
that do not contribute to the isospin asymmetry. Furthermore we discard corrections 
to 6| that are sub-leading in Ky^jmB- The calculation of such effects requires a better 
understanding of the light-cone wave functions of the B meson and is beyond the scope of 
this paper. However, we have convinced ourselves that the analogous effects that enhance 
the isospin asymmetry for transverse polarizations are absent in the longitudinal case, 
see discussion after (|14D below. 



2.2.1 Details of the calculation 

(a) Annihilation graphs: The annihilation graphs in Fig. ^(a) involve hadronic ma- 
trix elements of the four-quark operators Oi_6- To first approximation, these can be 
factorized into two matrix elements defining meson decay constants and meson-photon 
transition form factors, according to the two graphs in Fig. 0(a). The meson-photon 
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transition form factors and their role in the annihilation contribution to B decays have 
been extensively discussed in the recent literature |^-|^ . 



Let us start with the matrix element of the operators and Oq which can be 
expressed in terms of the (pseudo-) scalar K*'j* transition form factor and the B meson 
decay constant Jb- (An analogous term with the (pseudo-) scalar B'j* transition form 
factor is absent, because the remaining matrix element of a pseudo-scalar current between 
a vector meson and the vacuum vanishes.) Note that the annihilation graphs with O^^q 
are independent of the B meson wave function at the considered level of approximation. 
We calculate the (pseudo-) scalar K*'j* transition form factor perturbatively in terms of 
the K* light-cone wave functions, which yields 



{K*{p\eh*{q,fiMl + ^,)q\0) 

u m\ + uq^ 

+ {Cq ^ Cs, U^U) . (14) 



2 ^ e e, /i, ((g ■ e*) p'^ - e*^ (q ■ p') - le^^^^ eUpp'^) £ du ■ 



Note that only transverse polarizations do contribute here. (This remains true if one 
includes sub-leading contributions from twist-3 wave functions. Therefore the part of the 
isospin asymmetry which is related to longitudinal K* polarizations receives no enhanced 
contributions from the Wilson coefficients C5 and Cq.) Comparing with (|ll|) we obtain 
the following contribution to the function 

Ki^''\<f) = - + ^] F^iqVml) , F^{s) = \ f du ^ , (15) 

\ I\C J 6 Jo U + us 

which coincides with the result in for —* 0. (When comparing with |^ one 
has to use Ac ~ — Aj. Note that in ([l5|) and in the following we keep terms with s 
in the denominator of convolution integrals with light meson distribution amplitudes, 
although these terms are formally suppressed by s Ins for small values of s = g^/m^.) 
The functions K2 and do not receive contributions from C5 and Cq to the considered 
order in the A/j/m^ expansion. 

The {V—A)®(y—A) currents from the operators Oi-O^ give the leading contribution 
to the function k\ quoted in Eq. (0). For transverse K* polarizations one obtains the 
contribution 

K^'''\<f ) = (t +^^) + (^^ + f ) ' 

which again coincides with the term quoted in |2^. Note that the separate terms in 
the decomposition of {K*'~f*\{V — A) ® {V — A)\B) are invariant under QED gauge 
transformations only up to a "contact term" which does not depend on the meson's 
wave functions [|38|, In the sum, of course, these terms cancel. 

We finally mention that there is another mechanism to induce the B — >■ K*'~f* decay 
where the photon does not couple in a point-like way, but enters through its hadronic 
component [|5, I^. This effect gives rise to l/mj, corrections to the meson-photon 
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transition form factors, which are not taken into account here. We will also neglect NLO 
order terms arising from corrections to the annihilation topologies. In the numerical 
analysis below, we will use leading-order values for the Wilson coefficients Ci and C2, 
and next-to- leading order values for C3 — Ce, and include the uncertainty with respect 
to variations of the factorization scale /i into our error estimate. 



(b) Hard spectator interactions involving Og: Calculating the decay amplitude 
from the diagrams in Fig. @(b), with the insertion of the gluonic penguin operator O^, 
we obtain 



FAs) 



1 

3 Jo 



BJ 5 
1 



du ■ 



u + us)'- 



(17) 



(18) 



which again coincides with p6[ in the limit q 0. For the longitudinal part we have 



B) 5 



du 



011 (^) 
U + US 

'(1), 



(19) 



which has already been calculated in ||22|. Note that the function X]_ [s] contains a 
logarithmic divergence for s ^ which indicates the breakdown of QCD factorization 
for matrix elements with at sub-leading order in the A/j/m^ expansion. We follow 
the authors of |^ and parameterize the divergent part of the integral by replacing 



du 



pe 



i(f>\ 



l-Aft/ms 



du 



(20) 



and taking Ah — 0.5 GeV as a regulating cut-off, together with < p < 1 and 0, giving 
a conservative estimate of the theoretical uncertainty related to the absolute value and 
phase of this non-factorizing soft contribution. 



(c) Hard spectator interactions involving Oi^q. The calculation of the diagrams 
in Fig. ^(c) splits into two parts: First one calculates the internal quark loop which 
results in an effective b ^ sg* vertex with vector- and tensor-like couplings and form 
factors depending on the internal quark kinematics and the quark masses in the loop. 
The tensor contribution is local, and can be, as usual, absorbed into Cf^. The vector 
form factor is expressed in terms of the function [z = Am'^/s) 



his, m,) = -- (^In -f - - - .j - - (2 + .) 



1 

arctan , z > 1 

1 + z in 

m ^ z < 1 

x/i 2 - 

(21) 
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and reads 

Fv{s) = - {h{s, (C2 + + Ce) + h{x, m^) (C3 + C4 + Ce) 

+h{s, 0) (C3 + 3^4 + 3Ce) - ^ (C3 - - ISCe)} . (22) 

The rest of the calculation is completely analogous as for the diagrams with Og in 
Fig. 0(b). This yields the following contributions to the isospin asymmetry, 

Kt'-Hl') = ^tI t'i-^T^,Fv{u< + uf) : (23) 

A^c 47r 3 JO u + uq^ 

K-^^\q^) = -%^\ [ du (.) - ^1^) Fy[uml + uq^); (24) 

for the transverse part, and 

K^i^^\q') = -^^2 £duMu)Fy{uml + uq') ; (25) 

for the longitudinal one. In the limit ^ we reproduce the result in p6[ (neglecting 
the small Wilson coefficients C3-Q in Fy(s)). 



2.2.2 Numerical estimate 



We now turn to the numerical analysis of the isospin asjTiimetry in B 
combining the above results as 



Ki{q') 
Ki{q') 
Kl{q') 



Ki^-\q^) + Ki'^'\q^) + Ki^'\q^) 
K^^^\q') + K^^^\q') , 
Kl^^\q') + KY^{q^) + Kl^^\q') . 



(26) 



The central values and uncertainties of all hadronic and SM input parameters are taken as 
in Ref . p2| . The renormalization scale in the annihilation and hard-scattering processes 
will be taken as /i' = y/ fiA^, where /i will be varied as mfc/2 < yU < 2mb- 

The result for the isospin asymmetry is plotted as a function of in Fig. ^ At 



(27) 



g = we basically recover the result of Kagan/Neubert p6| with 



Re 



where the uncertainty due to the theoretical input parameters (taken as in p2|]) is dom- 
inated by the soft form factor (~ -1,5%), the B meson decay constant (~ 1.5%), 
the moment of the B meson light-cone wave function (~ 1.1%), and the IR-sensitive 
integral X^^ (~ 1%). We thus confirm the conclusion in that QCD factorization 



correctly reproduces the sign and magnitude of the experimentally measured isospin 
asymmetry {Aj[B K*'j] = 0.11 ± 0.07 [0, ||, H) with SM values for the Wilson 
coefficients. 
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Figure 3: The SM isospin asymmetry dAj/dq^ for the decay B K*i^i as a function 
of g^. The error band includes the variation of hadronic and SM input parameters (taken 



as in p^]) and of the renormahzation scale /i between mb/2 and 2mb 



For increasing values of q"^ the isospin-asymmetry decreases, and its central value 
becomes slightly negative above = 2 GeV^ and stays basically at a constant value of 
about —1%. Since the uncertainty related to the hadronic input parameters is reduced 
as well, this means that the measurement of a significant deviation from zero of the 
isospin asymmetry in the range 2 GeV^ < < 7 GeV^ may still indicate new physics 
(although one would need to have a handle on even higher order effects, before drawing 
any definite conclusions). Note that in the SM the isospin asymmetry is sensitive to C5 
and Cq at small but to C3 and C4 at larger g^. Thus, in principle, the two momentum 
regions provide complementary tests of the four-quark penguin-operators (see also Fig. || 
below) . 



3 Decay asymmetries in MSSM with MFV 

In this section we present an example to illustrate how new physics may affect the FB 
and isospin asymmetry in 5 — ^ K*i~^i~ decay. We consider the minimal supersymmetric 
standard model (MSSM) with minimal flavor violation (MFV) |Ty, |^|. In this model, 
all flavour transitions are governed by the CKM matrix. In addition, in order to increase 
the predictivity of the model, we impose the following restrictions. 

• All supersymmetric particles are taken as heavy (about 1 TeV), except for charginos, 
sneutrinos, the light (mostly right-handed) stop, and charged Higgs fields. Heavy 
particles (squarks, gluinos etc.) are integrated out to obtain a "low-energy" effec- 



tive theory in terms of light SUSY particles and SM fields [10 



The weak effective Hamiltonian (|I]) includes only SM operators. 

We assume a flavor diagonal mass matrix in the down-squark sector |p . Therefore 
neutralino contributions to the Wilson coefficients will be absent. 
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We carefully consider the two regimes of small and large tan /?, taking into account the 
modifications and subtleties associated with large values of tan (3 Jll], |5^ (tan j3 is the 
ratio of vacuum expectation values of the two Higgs doublets). 

The analysis is organized in three steps: First, we perform a scan of model parameters 
within specific ranges which will be indicated below. For each set of parameter values 
we will calculate the Wilson coefficients contributing to the B —>■ Xsj decay at NLO 



following |T0| , |TT], including also electroweak corrections and non-perturbative 
contributions [^J. Then we will use the experimental data on B —>■ Xg'-f to constrain 
the model parameter space, i.e. those parameter sets that lead to Wilson coefficients 
which are incompatible with the B Xg'j constraint will be discarded. (Note that the 
allowed region of parameter space will divide into two classes, corresponding to negative 
or positive sign of C^^.) Finally, using the remaining part of parameter space we evaluate 
the Wilson coefficients in (|I]), in terms of which we obtain numerical results for the two 
decay asymmetries, dAp^ and dAj in MFV. 

3.1 Model set-up 

The free parameters of the model are: 

• The Higgs mixing parameter /i. We consider |/i| < 500 GeV. 

• The mass parameter M2 in the chargino mass matrix ("weak gaugino mass"), 
which is varied as M2 = 100 — 800 GeV. Alternatively, one could consider one 
of the chargino mass eigenstates. In the parameter scan we always impose the 
constraint that the masses of both charginos should be heavier than 100 GeV, as 
required by direct searches. 

• The light stop mass is taken in the range 100 < rht2 < 350 GeV. (Note that 
depending on the relative size of stop, gluino and other squark masses, one has to 
be careful when evaluating the Wilson coefficients. In our case, we are considering 
a scenario with a very light right-handed stop and charginos in comparison with 
gluino and other squark masses. This implies that we have to use the Wilson 



coefficients at the weak scale as described in Eq.(31) of [|Tl|]. 



\9t\ is the stop mixing angle, such that ti = cos6ttL + sinOttR and ^2 = sinOttL + 
cosOftfi. Constraints from electroweak precision measurements favor the light stop 



quark to be almost right-handed We therefore take \6t\ < tt/IO. 



The charged Higgs mass is considered in the range Mh = 100 — 400 GeV. We 
have chosen a not too heavy mass for the charged Higgs to evaluate the maximal 
possible impact in our observables of a non-decoupled charged Higgs. 

The ratio of vacuum expectation values for the two Higgs doublets is varied as 
tan/3 = 2 — 40. In the expressions for the Wilson Coefficients the large tan/3 



regime is treated as in Ref. [|Tl|]. Large values of tan/? are particularly interesting 

'7 



since this can easily lead to a sign flip in the value of Cf^ 



The following parameters will be kept fixed in the analysis: 
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According to the effective theory approximation of the MFV supersymmetric model, 
all heavy squark masses (including the heavy stop) are taken at a common value, 
mq=l TeV. Also the gluino mass is taken at a fixed value iTLg=l TeV. The renor- 
malization effects from integrating out squarks and gluinos between the scale 1 TeV 



and the electroweak scale as discussed in 11 are taken into account 



- For the sneutrino masses we consider a conservative value, mj>=100 GeV which is 
not in conflict with direct searches. 

Scanning over discrete points in the so-defined parameters space, we calculate the Wilson 
coefficients for the operators in (|l|) at the electroweak scale, = 0{Mw). At NLO 
accuracy one obtains the generic form 

= C'f\l^w)sM + ^Cf\l^w)H + ^Cf\lJ,w)sVSY 



+ [Ci MsM + SCl^' Mh + ^CT' Msvsy) , (28) 

where we indicated explicitely the contributions from different virtual particles in box 
and penguin diagrams: 

• SM contributions from charged gauge bosons and up-type quarks. 

• Contributions from charged Higgs bosons and up-type quarks, present in SUSY, 
but also in generic models with two Higgs doublets. 

• Specific SUSY contributions from charginos and up-type squarks, gluinos and 
down-type squarks, neutralinos and down-type quarks. (We will mainly focus 
on chargino effects. In our model gluinos are very heavy and their main effect is 
to induce renormalization effects in the vertices involving charginos, charged Higgs 



bosons, charged gauge bosons and unphysical scalars |T^. As stated above, neu- 
tralino contributions are absent due to our assumption about the mass matrix in 
the down-quark sector.) 



The expressions for 6C^%^^'^^ can be found in Ref [0; SUSY contributions to Cg and 
Cio are taken from [1^, [1^. We have also included the NLO SUSY contribution to the 
four-quark penguin operators from Ref. . (Notice the different conventions used for 6t 



and the sign of the fi parameter in [|10| and 0.) The values of the Wilson coefficients are 
evolved to the low scales relevant for the exclusive B decays by standard renormalization 
group analysis. 

3.2 Constraints from B Xg^ 

We consider the world average of measurements of the inclusive branching ratio for 



2.82 X 10"^ < BR{B X,7) < 3.62 x 10"^ . (29) 
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The precision reached in the determination of this decay mode significantly constrains 
new physics contributions and requires to consider very carefully higher order QCD 
effects in the theoretical calculations, also in the presence of new physics. For instance, at 
two loops the electroweak amplitude gets enhanced in the SM due to a large QCD mixing 
with the effective operator b — >■ sec. Moreover, in the MFV model we are discussing it 
is known that there is a compensatory effect between the charged Higgs and chargino- 
stop contribution that gets unbalanced once the NLO correction is included. Apart 
from the standard renormalization group improvement of Wilson coefficients at NLO, 
we finally have to include non-perturbative corrections associated with the hadronic 



matrix elements for B Xg'j. In the analysis we will strictly follow Refs. [jl^, • 

Scanning over the allowed parameter space (notice that in our scenario we are only 
exploring light stop, charginos and not too heavy charged Higgses), we find the following 
allowed range for the new contributions 6Cj to the Wilson coefficient C^^ at the scale 
^iw (see Eq. (g8|)) 

- 0.02 < SCrifXw) < 0.12 or 0.59 < SCjifXw) < 1-24, (30) 

where the first range is SM-like (negative sign for C^^), and the second range corresponds 
to the fiipped-sign solution. The behaviour of the two ranges as a function of the param- 
eters tan/5, 6t and is quite different. We have found that the points in the SM-like 
region corresponds predominantly to values for tan f5 in the low regime (between 2 to 10). 
Moreover, higher values of rrih tend to prefer lower values of tan (3. There is no preferred 
value for 9t and also any value of Tht2 is found. The flipped-sign solutions for Cj^ corre- 
spond to large values of tan/5 of order 40. In this case, 6t preferably takes extreme values 
around ±7r/10. The allowed values for mt2 increase with increasing ruh but are always 
below 300 GeV, for example, for m/j = 100 GeV one gets 171^2 = 100 — 240 GeV, and 
for the highest value ruh = 400 GeV one has Tht2 = 100 — 300 GeV. The corresponding 
range for the Wilson coefficient Cf^ reads 

-0.18<5C8if^w) <0m or 0.20 < SCsifXw) < 0-^7 , (31) 

Concerning the Wilson coefficients C3 — Cq of four-quark penguin operators our MFV 
model is SM-like. Also the Wilson coefficients Cg and Cio only receive tiny contributions 
(of at most 10% from their SM value, we restricted ourselves to LO SUSY contributions 
here). The mechanisms that enhance the SUSY contributions to C^^ at large tan/? are 
not working for Cg and Ciq. For instance, the charged Higgs contribution dominant for 



Cj at large tan/? is suppressed for Cg^io fl^, Therefore the modifications for the FB 



and isospin-asymmetry, to be discussed below, will be due to new physics contributions 
to Cf, Cf mainly. 

3.3 FB and isospin asymmetry in MFV 

Using the parameter values in the MFV scenario allowed from the constraints discussed 
in the previous subsection, we calculate the FB and isospin-asymmetry as explained 
in section 0. The result for the FB asymmetry is shown in Fig. 11. The two distinct 



bands refer to the two possible signs of Cj in (|30|) . Comparing with the SM case in 



13 



-0.2 
-0.3 



C7>0 
C7<0 



3 4 

«2 



Figure 4: FB asymmetry (§) for the decay B K*i^i~ as a function of g^. The hght 
band corresponds to the region of parameter space with same (negative) sign for C7 as 
in the SM; the dark band refers to solutions with flipped sign for C7; and the dotted line 
is the central value of the SM prediction. 



Fig. |I], we observe that the variation of MFV parameters for the case of SM-like sign 
of C^^ (light band) is indistinguishable from the hadronic uncertainties. An exception 
is the FB asymmetry zero where the hadronic uncertainty is reduced due to the form 



factor relations discussed in p3|, ^ Here the MFV effects are competitive with the 
hadronic uncertainty, and show a tendency towards lower values of the FB asymmetry 
zero. This implies that if the experimental measurement of the FB asymmetry zero found 
an even larger values than predicted in the SM it would not easily be accommodated by 
our MFV scenario, and eventually would lead to the exclusion of the region of parameter 
space described in section [Ol . As is well known, for the flipped sign solution of C^^ (dark 
band) the characteristic features of the FB asymmetry change dramatically: dAps/dq^ is 
positive for small and large values of g^, and there is no asymmetry zero in the spectrum. 

The MFV prediction for the isospin asymmetry, shown in Fig. ^, also splits into two 
bands, according to the sign of C^^ . Because the isospin asymmetry is dominated by the 
penguin coefficients Cs-e which are only slightly modified in our MFV scenario, both 
bands are very thin. Within the MFV the isospin asymmetry at small values of thus 
provides a possibility to exclude the fiipped sign solution for Cf^ , independently from 
the analysis of the FB asymmetry (as already mentioned, at present, experimental data 
favor the same-sign solution). On the other hand, a measurement of a large isospin 
asymmetry at moderate values of would point to physics beyond both, the SM and 
MFV. 



4 Beyond MFV 

As we have seen in the previous section, in MFV the main effect on the FB and isospin 
asymmetry is due to a possible sign flip in the value of C^^ . Taking the central value of 
the experimental measurement of the isospin asymmetry in i? ^ K*'-) at face value, one 



14 




Figure 5: Isospin asymmetry (j^) for the decay B K*i~^i as a function of g^. Notations 
as in Fig. ^. 




Figure 6: The isospin asymmetry dAj/dq"^ for the decay B K*£^i as a function of g^. 
(a) The combination ag^"* = {Cq + 6*5/8) in the function K^^""^ is varied within a factor of 
two around its SM value (dark band = larger values, light band = smaller values), (b) 
The same for the combination a^^^ = {C4 + C^/3) in the functions kI^""^ and K2^"'\ 



is tempted to consider this scenario as already excluded. (A reduction of theoretical and 
experimental uncertainties for this observable, combined with a measurement of the FB 
asymmetry would allow for an ultimate conclusion.) The Wilson coefficients for semi- 
leptonic and 4-quark penguin operators are only slightly affected in our MFV scenario, 
and therefore for SM-like sign of C^^ MFV is indistinguishable from the SM within the 
theoretical uncertainties. 

On the other hand, in a generic new physics model, one might expect sizeable correc- 
tions to other Wilson coefficients as well, which may show up in the B — »• K*£~^£~ decay 
asymmetries. In the following we perform a brief discussion on the possible size of such 
effects. 

• The penguin operators Os-e give the main effect to the isospin asymmetry in i? — >■ 
K*'y and B K*£^i~. More precisely, dAi/dq^ is sensitive to the combinations 
af = {Cq + 65/3) and af = {C4 + C3/3) entering the functions i^i^^"^ and /sTp^ 
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that describe the dominant annihilation contribution to the isospin asymmetry at 
small and large g^, respectively. In principle, the size of the corresponding Wilson 
coefficients is already constrained from many sources, for instance non-leptonic B 
decays, life-time differences of B mesons, kaon decays etc. Since the sensitivity 
to individual Wilson coefficients in these decays is not very high, and considering 
the possibility of some conspiracy between different new physics effects one may, 
however, still allow for sizeable deviations of Cs-e from their SM values. Fig. 
shows the effect on the isospin asymmetry if we allow for a variation of ag°'' and 04 
within a factor of two. In this C3iSG, clS already mentioned above, modifications of 
Og'^'* (respectively a^^^) are competitive with the theoretical uncertainties at small 
(respectively large) values of (at large values, however, the isospin asymmetry 
itself is rather small, and probably the necessary accuracy to see modifications of 
af^ will be hard to achieve.) 

The semi-leptonic operators 09,10 are relevant for B —>■ K*i~^i~ at not too small 
values of g^. As already stated above the FB asymmetry provides a rather clean 
test of the Wilson coefficients Cg^io (compared to this the isospin asymmetry is 
less sensitive to Cg and Cio, unless the experimental measurement leads to an 
unexpectedly large asymmetry for 2 GeV^ < g^ < 7 GeV^.) In particular, when 
NLO QCD effects are included the Wilson coefficient Cg can be measured with 10% 
accuracy from the FB asymmetry zero (provided the relative signs of 6*7,9,10 
are SM-like). If the relative sign between C^^ and Cg is flipped there will be no FB 
asymmetry zero at all; in case of a non-SM-like sign of Ciq the overall sign of the 
FB asymmetry will change. We refrain from giving quantitative results here, but 



rather refer to [Ig, |2J, ^ for extensive phenomenological analyses of these issues. 



In this work we do not entertain the possibility of introducing new operators in 
the effective Hamiltonian (|l]). An example of contributions of non-standard oper- 
ators to the FB asymmetry in i? ^ K*i'^i~ can be found in [Q. A discussion 
of non-standard penguin operators, which would also be relevant for the isospin 
asymmetry, has been given in |BB . 



5 Summary 

We discussed two dedicated observables in B ^ K*i'^i~ decays, the forward-backward 
and isospin-asymmetry, which can be used to test the standard model against new 
physics. We concentrated on the large recoil region where the decay amplitudes can 
be calculated using the QCD factorization approach. (For the FB asymmetry we used 
the results from for the isospin asymmetry we generalized the B — > K*'~f analysis 



performed in |^.) The two asymmetries turn out to give complementary tests of the 
Wilson coefficients in the weak effective Hamiltonian: For a given value (and sign) of 
C^^, the FB asymmetry is sensitive to the semi-leptonic coefficients Cg and Cio. On the 
other hand the isospin asymmetry is dominated by the effects of four-quark coefficients 
C5,6 and 6*3,4 at small or medium values of g^. 
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To illustrate possible new physics effects we considered an explicit model, namely 
minimal supersymmetry with minimal ffavor violation. In this scenario the main devi- 
ation from the SM is the possibility to generate another sign for the magnetic penguin 
coefficients Cj^, Cf^ (while the absolute value for C^^ is constrained by the B — > Xs'j 
measurement). The corresponding two regions of allowed parameter sets translate into 
two well separated bands of predictions for the FB and isospin-asymmetry. On the other 
hand, the change in the remaining Wilson coefficients is rather small in the MFV scenario 
and indistinguishable from hadronic uncertainties of the theoretical calculation. 

Considering an alternative point of view and assuming the SM to be a valid approx- 
imation in rare B decays, the simultaneous study of the isospin asymmetry and the 
forward-backward asymmetry in the large- recoil region (g^ < 4m^) provides a perfect 
testing ground for the QCD factorization approach. 
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